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Definition 

Distance  of  pitch  axis  from  midchord  of  hydrofoil  model, 
in  units  of  semichord  (positive  aft) 

Local  aspect  ratio;  square  of  distance  from  local  position 
to  tip  of  hydrofoil  divided  by  one-half  the  area  from  local 
position  to  tip  of  hydrofoil 

Semichord  of  hydrofoil  model 

Bending  stiffness 

Torsional  stiffness 

Translation  of  foil  elastic  axis  from  unstressed  position 
at  station  i 

Reduced  frequency;  ba>/V 
Bending  stiffness  matrix 
Torsional  stiffness  matrix 
Length  of  hydrofoil  from  root  to  tip 

Steady  moment  on  foil  due  to  angle  of  attack  at  station  i 
Spanwise  modification  factor  for  noncirculatory  loading 
Steady  lift  on  foil  due  to  angle  of  attack  at  station  i 
Radius  of  gyration,  in  units  of  semichord 
Time 

Divergence  speed  of  foil 
Water  speed 

Distance  of  center  of  gravity  from  elastic  axis,  in  units 
of  semichord  (positive  aft) 

Angle  of  attack  of  foil  relative  to  free  stream  at  station  i 

Angle  of  unstressed  foil  relative  to  free  stream  at  station  % 

Exponential  decay  factor  corresponding  to  oscillation  of 
amplitude  proportional  to  e~ 
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Definition 


Component  of  exponential  decay  factor  due  to  hydro- 
dynamic  damping 

Component  of  exponential  decay  factor  due  to  structural 
damping;  ^s^/u^T~sJ 

Structural  damping  parameter 

Torsional  displacement  of  foil  elastic  axis  from  unstressed 
position  at  station  i 

Mass  ratio:  ratio  of  mass  of  hydrofoil  to  mass  of  cylinder 
of  water  circumscribing  chord  of  foil 

Circular  frequency  of  oscillation 


ABSTRACT 


Four  low  mass  ratio  hydrofoil  models  of  aspect  ratio  5  were  flutter 
tested.  The  flutter  speed  of  the  mass  ratio  0.963  model  was  24.7  knots. 

The  other  three  models,  of  mass  ratios  0.455,  0.395,  and  0.202,  were  sub¬ 
ject  to  static  failure  in  bending  at  about  36  knots  and  did  not  flutter  below 
this  speed.  The  results  were  compared  with  the  predictions  of  three  flutter 
theories.  All  theories  gave  conservative  flutter  speed  predictions  at  mass 
ratio  0.963.  Two  of  the  theories  were  also  conservative  at  mass  ratio  0.455. 

The  other  predictions  could  not  be  evaluated.  Model  divergence  character¬ 
istics  were  responsible  for  model  failure  in  bending. 

ADMINISTRATIVE  INFORMATION 

The  work  performed  at  the  Naval  Ship  Research  and  Development  Center  was  funded  by 
the  Naval  Ship  Systems  Command  under  Subproject  S4606,  Task  1703. 

The  work  performed  at  Southwest  Research  Institute  was  sponsored  by  the  Hydro¬ 
mechanics  Laboratory  of  the  Naval  Ship  Research  and  Development  Center  under  the  Naval 
Ship  Systems  Command  Subproject  S4606,  Task  1703.  Funding  was  provided  according  to 
Contract  N00014-69-C-0219. 


INTRODUCTION 

A  long-standing  difficr’ty  in  the  field  of  hydrofoil  flutter  research  has  been  the  inabil¬ 
ity  of  flutter  theory  and  experiment  to  establish  the  dependence  of  hydrofoil  flutter  character¬ 
istics  on  mass  ratio  in  the  mass  ratio  range  below  1.0.  While  flutter  theory  has  generally 
failed  to  predict  experimental  results,1,  2'  3  flutter  experiments  as  well  have  often  failed  to 
produce  usable  results  in  the  form  of  flutter  occurrences,4, 5  In  cases  where  some  agreement 
has  been  obtained  between  theory  and  experiment,  the  data  have  been  insufficient  to  confirm 
theoretical  predictions  over  a  range  of  mass  ratio  values6  and  the  theories  have  lacked 
genera)  applicability3  or  self-consistency.1 

A  study  by  Rowe3  using  highly  accurate  structural  calculations  indicated  that  the  defi¬ 
ciencies  of  flutter  theory  lay  in  the  hydrodynamic  loading  formulation.  By  using  an  improved 
aerodynamic  strip  theory  in  a  six-mode  R.iyleigh-type  analysis,  Yates7  successfully  predicted 
the  flutter  characteristics  of  airfoils  of  mass  ratio  as  low  as  1.0  and  obtained  results  within 
20  percent  of  experiment  for  a  hydrofoil8  of  mass  ratio  0.99.  Yates  also  predicted  a  flutter 
boundary  at  very  low  mass  ratio,  characterized  by  a  flutter  speed  which  approached  zero  as 
mass  ratio  approached  zero.  The  existence  of  such  a  flutter  boundary  would  be  highly  signi  ¬ 
ficant  for  full-scale  hydrofoils. 
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1*  wa-  felt  that  Yates'  flutter  theory  was  a  significant  improvement  over  existing  flutter 
theories  in  its  modi fi eat  ion  of  two-dimensional  strip  theory  to  include  three-dimensional  ef¬ 
fect  -  nn  loading.  However,  additional  data  were  required  to  confirm  its  predictions  of  flutter 
at  low  mass  ratio.  Therefore  the  present  work  was  undertaken. 

Four  hydrofoil  flutter  models,  and  a  model  support  system,  were  designed  and  construct¬ 
ed  by  Southwest  Research  Institute  (SwRI). 9  The  flutter  models  ranged  in  mass  ratio  from  0.963 
to  0.202,  with  the  0.963  model  being  a  dynamically  scaled  model  of  the  mass  ratio  0.99  hydro¬ 
foil  mentioned  above.  The  models  were  flutter  tested  in  the  36-in.  variable  pressure  water 
'unnel10  at  the  Naval  Ship  Research  and  Development  Center  (NSRDC)  during  February  1970. 

The  flutter  characteristics  of  the  four  flutter  models  were  calculated  using  three  flutter 
theories.  F.ach  of  the  flutter  theories  incorporated  some  spanwise  modification  of  hydrodynamic 
loading,  hut  none  was  precisely  the  same  as  Yates’  original  calculation.  In  a  calculation 
made  at  NSRDC,  Yates’  loading  was  combined  with  a  lumped  parameter  structural  representa¬ 
tion  similar  to  Rowe.  Second,  the  ioad  theory  was  extended,  as  suggested  by  Yates,7  to  in¬ 
clude  three-dimensional  effects  on  noncirculatory  loading  and  was  used  with  the  lumped  pa¬ 
rameter  structural  representation.  The  third  calculation,  performed  atSwRI,  used  measured 
values  of  -panwise  lift  slope  in  a  two-mode  Rayleigh-type  analysis.'*  Results  calculated 
for  the  mass  ratio  0.99  hydrofoil  by  the  two  NSRDC  theories  are  also  presented.  Results  for 
this  hydrofoil  calculated  by  the  SwRI  theory  have  been  previously  published.11 

FLUTTER  MODELS 

Model  design  and  construction  have  been  fully  described  in  a  previous  report.9  The 
following  discussion  is  intended  as  a  summary  of  model  design  philosophy  and  characteristics. 

The  four  flutter  models  were  intended  to  form  a  family  in  which  mass  ratio  varied  to 
substantially  below  1.0  while  all  other  parameters  remained  constant.  Geometrical  size  and 
design  flutter  speeds  were  based  on  the  capabilities  of  the  test  facility,  the  36-in.  variable 
pressure  water  tunnel10  at  NSRDC.  A  parent  model  which  had  been  previously  flutter  tested8 
provided  values  for  the  invariant  structural  parameters  and  for  the  mass  ratio  of  the  heaviest 
of  the  new  models  (Model  1).  This  correlation  permitted  a  comparison  of  scaling  effects  be¬ 
tween  the  parent  model  and  the  new,  smaller  Model  1.  The  remaining  models,  Models  2,  3, 
and  4,  were  of  lower  mass  ratio  than  Model  1. 

The  parent  model  was  an  aspect  ratio  5  rectangular  hydrofoil  with  a  30-in.  span  (to  a 
reflecting  plate)  and  a  mass  ratio  of  0.99.  It  had  undergone  structural  failure  after  flutter 
occurred  at  48.1  knots,  with  a  frequency  of  17.5  Hz  and  a  reduced  irequency  of  0.676.  In 
order  to  ,  reserve  similarity  between  the  hydrodynamic  forces  acting  on  the  parent  and  Model  1, 
the  same  value  of  reduced  frequency  at  flutter  was  sought  for  Model  1.  The  50-knot  maximum 
speed  of  the  water  tunnel  and  its  36-in.  diameter  jet  led  to  a  choice  of  &ModeI  j  =  1/2  5parent 
an(i  *Modei  r  ^Parent’  which  would  hopefully  result  in  the  relationship  VModel  j=  1/2  Fparen( 
ut  flutter.  These  scale  ratios  gave  a  wide  speed  range  for  testing,  enabled  the  resulting 
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15-in.  models  and  a  reflecting  plate  to  be  located  well  inside  the  water  jet,  ar.d  represented 
feasible  structural  stiffnesses  with  El  and  GJ  being  1/64  of  the  parent  model  values.  The 
structural  parameter  values  of  the  parent  model  are  given  in  Table  1. 

The  models  were  constructed  by  cementing  airfcil-shaped  segments  to  steel  spars  at 
the  quarter  chord  position.  The  segments  were  ballasted  to  produce  the  desired  inertial 
characteristics.  Four  identical  spars  were  manufactured  from  high  strength  maraging  steel 
with  a  yield  point  of  350,000  psi.  Segments  were  individually  cast  of  the  following  materials 
to  produce  the  desired  variation  in  mass  ratio: 


Model 

Segment 

Material 

1 

73  percent  lead,  27  percent  tin 

2 

epoxy  and  lead  powder 

3 

epoxy  and  lead  powder 

4 

epoxy  and  glass  microballoons 

One  of  the  models  is  shown  in  Figure  1  prior  to  assembly.  Each  model  contained  an  exciter 
rod  anchored  to  its  spar  tip,  passing  through  the  center  of  the  spar,  and  extending  beyond  the 
root  end  of  the  spar.  Strain  gages  sensitive  to  bending  and  torsion  were  cemented  to  the  root 
of  the  spar  in  a  hollowed-out  portion  of  the  root  segment.  After  assembly,  gaps  between  seg¬ 
ments  were  sealed  with  Silastic  RTV  731  compound  to  form  a  smooth  surface.  Two  of  the 
completed  models  are  shown  in  Figure  2. 

Inertial  parameter  values  for  the  completed  models  are  given  in  Table  1,  along  with  the 
values  for  the  parent  model.  Values  of  mass  unbalance  xa  and  radius  of  gyration  rfl  agreed 
extremely  well  with  design  values.  The  mass  ratio  of  Model  1  was  slightly  below  the  intended 
value  of  0.99,  while  the  mass  ratio  of  Model  4  was  the  lowest  obtainable  with  available 
materials. 


MODEL  SUPPORT  SYSTEM 

Each  model  was  supported  for  flutter  testing  by  an  enclosure  which  was  bolted  to  the 
ceiling  of  the  water  tunnel,  as  shown  in  Figure  3.  The  bottom  of  the  enclosure  consisted  of  a 
reflection  plate  40  in.  long  and  24  in.  wide  the  upper  surface  of  which  was  bevelled.  The  re¬ 
flection  plate  extended  3  in.  into  the  water  jet,  and  the  model  extended  15  in.  further  to  the 
centerline  of  the  water  jet.  The  enclosure  had  an  NACA  16-022  profile. 

Within  the  enclosure  were  housed  mechanisms  for  adjusting  the  angle  of  attack,  or  trim, 
of  the  flutter  model  and  for  twisting  the  exciter  rod  in  the  model  to  produce  a  tip  deflection. 

The  model  was  clamped  in  a  turntable  mounted  flush  with  the  bottom  of  the  reflection  plate,  as 
shown  in  Figure  4.  Trim  adjustment  was  made  by  rotating  the  entire  turntable  and  retainer 
plate  assembly  on  Teflon  coated  surfaces.  The  exciter  rod  was  twisted  by  a  cam  acting  on  the 
exciter  crank  independently  of  the  turntable,  allosving  the  tip  deflection  to  be  released  quickly 
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after  touching  its  maximum.  By  operating  the  cam  drive  motor  both  clockwise  and  counter¬ 
clockwise,  two  tip  deflection  amplitudes  were  obtained.  Initially  the  deflections  were  5  deg 
and  deg.  During  the  test,  the  earn  was  machined  down  to  give  deflections  of  2  deg  and 
2  .1  deg. 

This  method  of  deflecting  and  releasing  the  model  has  been  previously  used  to  measure 
damping  in  a  hydrofoil  flutter  test. 8  However,  the  hydrofoil  in  that  test  was  destroyed  when 
flutter  occurred.  It  was  felt  that  flutter  of  the  present  models  could  be  avoided  by  testing 
until  damping  had  become  very  small  and  then  removing  the  models  from  the  water  tunnel. 

The  flutter  inception  speed  would  be  determined  by  extrapolation  to  zero  damping. 

PROCEDURE 

STIFFNESS  MEASUREMENTS 

Bending  and  torsional  stiffness  measurements  were  made  prior  to  flutter  testing  on  all 
four  models  and  were  repeated  on  the  three  surviving  models  after  flutter  testing.  Pre-test 
measurements  were  made  with  each  model  clamped  in  a  vise,  while  post-test  measurements 
were  made  with  each  model  fixed  in  the  model  support  enclosure  which  was  used  during 
flutter  testing. 

The  models  were  statically  deflected  in  bending  or  torsion  in  the  following  manner. 
Bolts  were  placed  in  threaded  holes  in  the  tip  of  the  model,  one  in  the  center  of  the  spar  and 
the  other  2  in.  aft  of  the  spar  center.  Bending  deflections  were  produced  by  applying  loads 
at  the  spar  tip  perpendicular  to  the  foil  planform.  A  dial  gage  read  deflections  at  the  spar. 
Torsional  deflections  were  produced  by  placing  a  knife-edged  bolt  in  the  spar  tip  and  apply¬ 
ing  loads  to  the  other  bolt  while  the  knife  edge  rested  on  a  support.  A  maximum  of  15  lb  of 
weights  was  used  to  load  the  models. 

Bending  stiffness  El  was  calculated  from  the  relation12 


where  P  =  force  applied  to  tip  of  spar, 

y  =  spar  deflection  at  spanwise  position  x , 

L  =  length  of  spar,  and 

x  ~  spanwise  position  at  which  deflection  is  measured. 

This  relation  is  based  on  the  assumption  that  El  is  constant  along  the  spar.  The  torsional 
stiffness  GJ  was  calculated  from  the  relation13 
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where  M  =  moment  applied  to  tip  of  spar, 

x  =  spanwise  position  at  which  twist  angle  is  measured,  and 
d  -  twist  angle  at  spanwise  position  x. 

It  is  assumed  that  GJ  is  constant  along  the  spar. 


STRUCTURAL  MODE  MEASUREMENTS 


The  two  lowest  natural  frequencies  of  the  parent  model  were  determined  by  tip  excita- 
tior  at  the  time  of  its  flutter  testing.  The  frequencies  in  air  and  in  water  f.re  given  in 
Tab.'e  2. 

The  natural  frequencies  and  some  nodal  lines  for  the  current  models  were  obtained  by 
exciting  each  model  with  an  electromagnetic  shaker  placed  about  1/4  in.  from  the  surface  of 
the  model,  over  the  spar  and  near  the  model’s  tip,  as  shown  in  Figure  5.  There  was  no  me¬ 
chanical  connection  between  the  model  and  the  shaker.  The  shaker  was  driven  by  a  sweep 
oscillator  which  automatically  varied  its  frequency  over  a  preset  range.  Model  response  was 
registered  by  an  accelerometer  taped  to  the  model  on  the  opposite  side  from  the  shaker.  Ac¬ 
celeration  amplitude  was  automatically  plotted  as  a  function  of  frequency,  with  resonant  fre¬ 
quencies  corresponding  to  response  peaks.  A  second  “roving”  accelerometer  was  used  to 
determine  the  mode  shape  by  locating  nodal  lines.  Crossing  a  nodal  line  resuited  in  a  180- 
deg  phase  change  between  the  two  accelerometer  signals. 


Structural  damping  in  the  form  of  the  damping  constant  £ s  for  the  models  was  obtained 
in  two  ways.  In  the  first  method,  it  was  assumed  that  the  amplitude  of  decaying  oscillations 

~C  a>t 

was  proportional  to  c  5  .  The  constant  £  was  determined  directly  from  the  oscillations 

recorded  after  the  tip  of  the  model  was  twisted  and  released.  The  second  method  involved 


calculation  of  £s  from  the  width  of  the  model  response  curve  at  a  resonant  frequency. 
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FLUTTER  TESTING 

Flutter  testing  of  the  four  hydrofoil  models  was  performed  in  the  36-in.  variable  pres¬ 
sure  water  tunnel  at  NSRDC.  The  maximum  speed  through  the  test  section  is  50  knots. 

Static  pressure  at  the  test  section  centerline  can  be  varied  from  4  to  60  psia. 

The  model  support  enclosure  was  mounted  in  the  ceiling  of  the  open  jet  test  section, 
as  shown  in  Figure  3,  with  the  model  extending  to  the  centerline  of  the  test  section. 

After  installation  in  the  tunnel,  the  strain  gages  in  each  mode!  were  connected  to 
Endevco  Model  4402  signal  conditioning  units.  Gage  output  was  amplified  by  Dana  Model 
282J  amplifiers  and  recorded  on  a  Consolidated  Electrodynamics  Corporation  Model  5-124 
oscillograph.  The  model  was  excited  in  air  in  order  to  obtain  structural  damping  values  for 
the  two  lowest  modes. 

Since  both  the  enclosure  and  the  model  were  hollow,  it  was  decided  to  maintain  a 
slightly  higher  air  pressure  inside  them  than  in  the  test  section  to  prevent  water  from  entering. 


This  pressure  differentia!  had  to  bo  maintained  throughout  an  absolute  pressure  range  of 
about  9  to  -10  psia  which  occurred  during  filling  and  emptying  the  tunnel.  Following  the  test- 
sect  ion  pressure  at  all  times  would  also  prevent  possible  damage  due  to  pressure  loading  on 
the  enclosure  and  the  model.  A  pressure  regulation  system  was  constructed,  using  a  back¬ 
pressure  regulator  containing  a  teflon  diaphragm  as  the  controlling  element.  High  pressure 
air  was  bled  through  the  regulator  which  maintained  'he  desired  pressure  differential  by 
balancing  the  tunnel  pressure  against  the  enclosure  pressure. 

During  tunnel  operation,  air  constantly  leaked  out  of  the  enclosure  along  welded  seams 
which  had  not  been  completely  filled.  Air  leaks  also  occurred  from  the  Silastic  seams  in 
Flutter  Model  4,  because  lubricating  oil  from  the  exciter  rod  bearing  saturated  the  inside  of 
the  model  and  gradually  loosened  the  Silastic  bonding.  Several  repairs  were  made  to  the 
seals  in  order  to  complete  testing  Model  4.  Air  leaks  from  the  enclosure  resulted  in  a  con¬ 
tinually  increasing  air  content  in  the  tunnel,  primarily  evidenced  by  bubbles.  High  bubble 
content  might  affect  hydrodynamic  loading  in  addition  to  reducing  visibility.  Since  the  test 
section  was  emptied  each  day,  air  accumulated  for  a  maximum  of  seven  hours  while  testing, 
although  running  times  between  refillings  were  often  much  less  than  seven  hours.  No  esti¬ 
mate  of  the  effect  of  air  accumulation  could  be  made  from  the  test  data  obtained.  Some  re¬ 
duction  in  the  quantity  of  air  bubbles  was  achieved  by  increasing  the  static  pressure  in  the 
tunnel  from  20  psia  which  was  used  at  lower  speeds  to  25  psia  and  30  psia  for  higher  speed 
runs.  The  higher  pressures  also  reduced  a  small  amount  of  cavitation  which  occurred  along 
the  leading  edge  of  the  reflection  plate. 

Flutter  testing  was  performed  by  increasing  the  water  speed  from  zero  and  measuring 
the  model’s  damping  at  each  speed  selected.  After  zeroing  its  angle  of  attack,  c.,e  model  was 
excited  using  the  internal  exciting  mechanism.  The  exponential  decay  factor  S  was  then  cal¬ 
culated  from  the  recorded  signal.  Flutter  would  occur  at  S  =  0  so  that  a  decrease  in  S  toward 
zero  would  indicate  the  approach  of  a  flutter  instability.  This  method  successfully  detected 
the  approach  of  flutter  for  Model  1,  which  was  removed  from  the  tunnel  intact  after  coming 
within  an  estimated  1/2  knot  of  flutter.  It  was  found  helpful  in  the  low  damping  range  near 
flutter  to  make  several  damping  measurements  so  that  data  scatter  would  include  the  lowest 
damping  value. 

Of  the  two  excitation  amplitudes  available,  the  high-amplitude  excitation  was  used 
for  low  speeds  and  the  low-amplitude  for  high  speeds.  The  amplitude  was  reduced  when 
torsional  deflection  of  the  model  appeared  to  be  too  large.  Model  1  was  tested  at  5  deg 
excitation  up  to  17  knots,  and  at  2  deg  excitation  up  to  its  maximum  speed  of  24,2  knots. 
Models  2  and  4  required  a  reduction  in  amplitude  to  2/3  of  a  degree  for  speeds  above  30 
knots.  The  reduced  amplitude  was  obtained  by  machining  down  the  exciter  cam.  It  was 
found,  however,  that  the  2/3  deg  excitation  did  not  give  a  sufficiently  large  torsional  deflec¬ 
tion  for  the  damping  to  be  found.  This  difficulty  appeared  to  be  related  to  the  divergence 
characteristics  of  the  models  as  will  be  discussed  elsewhere. 
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Several  excitation  responses  of  the  models  were  recorded  on  high-speed  motion  picture 
film,  taken  at  500  frames/sec.  The  camera  photographed  the  motion  of  the  tip  of  the  model 
from  below  the  tunnel  test  section. 


RESULTS 


STIFFNESSES 

Bending  and  torsional  stiffness  values  are  given  in  Table  1,  along  with  the  other 
structural  parameters.  Values  measured  after  the  fluiter  test  agreed  with  those  measured  be¬ 
fore  the  test.  The  values  given  correspond  to  model  loading  in  the  high  load  range,  around 
i5  lb,  rather  than  the  low  load  range.  Higher  loading  gave  lower  stiffnesses  because  the 
structural  response  was  nonlinear.  The  hysteresis  curve  of  the  model  for  cyclic  loading  was 
approximately  represented  by  the  chosen  stiffnesses. 

The  design  values  for  El  and  GJ  were  53,100  lb-in.2  and  15,200  lb-in.2,  respectively. 
Therefore  Models  2,  3,  and  4  were  in  good  agreement  with  the  design  values,  while  Model  1 
exhibited  a  torsional  stiffness  which  was  somewhat  high. 

STRUCTURAL  MODES 

The  hydrofoils  exhibited  structural  vibration  modes  which  contained  both  bending  and 
torsion  components.  However,  the  two  lowest  modes  were  observed  to  be  primarily  first  bend¬ 
ing  and  first  torsion  modes,  respectively. 

Measured  natural  frequencies  are  given  in  Tables  3  through  6.  All  frequencies  up  to 
approximately  200  Hz  were  determined  by  forced  oscillation,  using  an  electromagnetic  shaker, 
but  only  the  first  two  modes  could  be  excited  by  tip  excitation.  Model  3  was  not  installed  in 
the  tunnel  and  therefore  no  tip  excitation  frequencies  are  available.  The  strain  gage  designed 
for  bending  on  Model  4  contained  a  broken  lead  wire  so  that  no  bending  frequencies  could  be 
obtained.  The  frequencies  obtained  by  the  two  methods  of  excitation  agreed  wi'hin  the  experi¬ 
mental  uncertainty  of  four  percent  for  mode  1,  but  did  not  agree  for  mode  2  of  Models  1  and  2. 

A  comparison  of  the  natural  frequencies  for  Model  1  and  the  parent  model  showed  that 
Model  1  had  frequencies  which  were  15  percent  and  7  percent  higher  than  the  parent  model  in 
modes  1  and  2,  respectively.  The  higher  frequencies  were  due  to  the  relatively  higher  stiff¬ 
nesses  of  Model  1,  given  in  Table  1,  which  exceeded  the  stiffness  scaling  ratio  of  1/64. 
Therefore  Model  1  was  not  a  perfectly  scaled  version  of  the  parent  model,  and  was  not  expect- 
ed  to  have  the  design  flutter  characteristics  of  a  24.05  knot  flutter  speed  and  a  17.5  Hz  flutter 
frequency.  However,  the  model  was  a  close  approximation  to  the  desired  scaling. 

Nodal  lines  obtained  during  forced  oscillation  are  show'n  by  the  data  points  in  Figure  6. 
The  points  represent  a  180-deg  phase  change  in  acceleration.  Theoretical  noda!  lines  are  also 
shown  in  Figure  6  and  will  be  discussed  later.  The  nodal  lines  for  Mode  1,  primarily  first 
bending,  could  not  be  detected  along  the  roots  of  the  models.  Modes  2  and  3  were  measured 
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on  each  model  except  Model  3,  for  which  the  nodal  line  of  mode  3  "ould  not  be  detected.  The 
nodal  lines  were  generally  smooth  curves  showing  only  slight  changes  from  model  to  model. 
Mode  2  of  Model  4,  however,  contained  an  irregularity  between  the  6  th  and  7  th  segments  from 
the  tip. 

Structural  damping  values,  obtained  from  tip  excitation  in  air  and  in  water  and  from  the 
forced  oscillation  response  curves,14  are  given  in  Table  7.  The  torsional  damping  of  Model 
4  in  air  could  not  be  determined  because  the  torsion  gage  output  contained  toe  large  a  bending 
component. 

The  two  methods  of  measuring  damping  did  not  agree,  although  damping  values  obtained 
by  either  method  were  fairly  consistent.  An  evaluation  of  the  experimental  uncertainty  in  the 
two  methods  indicated  that  there  was  considerably  more  uncertainty  in  the  damping  obtained 
from  forced  oscillation  response  peaks.  The  uncertainty  originated  in  plotting  and  measuring 
the  response  peak  heights  and  widths,  and  in  possible  variation  in  the  driving  force  of  the 
oscillator  as  the  model  response  varied.  Since  damping  values  obtained  from  twisting  and  re¬ 
leasing  the  models  were  considered  more  accurate,  these  values,  taken  from  in-water  model 
excitation,  were  used  in  the  flutter  calculations.  In  view  of  the  agreement  in  measured  in-air 
damping  between  Models  2  and  3,  the  in-water  damping  of  Model  2  was  used  in  flutter  calcula¬ 
tions  for  Model  3. 

FLUTTER  CHARACTERISTICS 

Models  1,  2,  and  4  were  tested  to  determine  their  flutter  characteristics.  Model  1  was 
tested  until  flutter  was  imminent  and  was  then  removed  for  the  tunnel.  Model  2  was  tested 
until  a  static  failure  in  bending  occurred,  with  no  evidence  of  flutter.  The  damaged  model  is 
shown  in  Figure  7.  Model  4  was  tested  and  removed  before  a  static  failure  occurred,  without 
evidence  of  flutter.  Because  of  the  similarity  between  Models  2  and  3,  Model  3  was  not  tested 
after  the  static  failure  of  Modei  2. 

Experimental  values  of  the  exponential  decay  factor  and  oscillation  frequency  for 
Model  1  are  plotted  in  Figures  8  and  9,  respectively.  The  highest  speed  attained  was  24.2 
knots.  In  view  of  the  rapidly  decreasing  decay  factors  it  was  estimated  that  flutter  would 
have  occurred  within  one-half  knot  of  the  highest  speed  tested,  so  that  a  flutter  speed  of 
24.7  knots  is  projected,  with  an  uncertainty  of  ^0.5  knot.  The  flutter  frequency  is  projected 
to  be  20.5  Hz,  which  was  the  lowest  frequency  that  occurred  during  testing  and  which  corre¬ 
sponded  to  the  lowest  decay  factor.  The  reduced  frequency  at  flutter  is  therefore  0.772, 
based  on  the  above  projected  values. 

The  design  flutter  speed  for  Model  1  was  24.05  knots,  with  a  flutter  frequency  of  17.5 
Hz.  Therefore  the  actual  flutter  speed  was  3  percent  higher  than  design,  while  the  flutter  fre¬ 
quency  was  17  percent  higher  than  the  design  value.  The  scaling  of  the  parent  model  was 
successful  in  that  flutter  speed  was  reduced  to  a  value  close  to  half  its  original  value,  while 
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the  discrepancy  in  flutter  frequency  can  be  explained  by  the  deviation  in  the  measured  stiff¬ 
nesses  and,  to  a  lesser  degree,  a  small  discrepancy  in  mass  ratio. 

The  results  for  Model  1  demonstrate  that  it  is  possible  to  approach  a  hydrofoil  flutter 
condition  very  closely  without  destroying  the  model.  Very  small  speed  increments  were  used 
to  obtain  a  total  damping  value  which  was  54  percent  of  the  structural  damping  or  1.3  percent 
of  critical  damping  when  testing  of  Model  1  was  stopped.  The  decrease  in  total  damping  to  a 
value  smaller  than  the  structural  damping  indicated  that  the  hydrodynamic  damping  had  al¬ 
ready  become  negative.  It  is  concluded  that  the  flutter  inception  speed  of  a  hydrofoil  can  be 
determined  in  some  cases  without  the  occurrence  of  flutter  if  damping  measurements  are 
made  while  speed  is  increased  and  if  speed  increments  are  kept  small. 

Model  2  was  tested  in  the  same  manner  as  Model  1.  However,  Model  2  displayed  ex¬ 
cessively  high  static  deflections  during  tip  excitation  in  the  speed  range  above  25  knots,  and 
eventually  failed  in  bending  at  36  knots.  Exponential  decay  factors  and  frequencies  were 
measured  using  a  maximum  of  2  deg  of  tip  twist  up  to  30  knots  and  are  shown  in  Figures  10 
and  11.  The  decay  factors  may  have  begun  to  decrease  at  30  knots,  but  no  further  measure¬ 
ments  were  possible  at  this  excitation  amplitude  because  of  the  alarming  deflections  of  tho 
model.  The  deflections  had  a  mean  value  which  gradually  decreased  to  the  unstressed  posi¬ 
tion  after  the  excitation  was  released.  During  this  decrease  a  few  cycles  of  oscillation 
occurred.  By  reducing  the  excitation  angle  to  2/3  deg,  the  model  was  tested  to  35  knots. 

The  smaller  excitation  produced  large  decay  factors  which  were  beyond  the  scale  of  Figure 
10  at  20  knots,  and  which  could  not  be  read  at  25  knots  and  above.  The  large  decay  factors 
were  caused  by  an  excitation  which  did  not  sufficiently  exceed  the  ambient  flow  fluctuations 
and  by  the  tendency  of  the  model  to  respond  excessively  at  higher  speeds.  Model  response 
to  flow  fluctuations  at  35  knots  was  visible  both  on  the  oscillograph  recordings  and  to  the 
naked  eye.  Since  damping  could  not  be  measured  under  these  conditions,  it  was  decided  to 
increase  tunnel  speed  gradually  while  recording  model  response  until  the  model  failed;  the 
failure  would  have  to  be  interpreted  as  either  flutter  or  static  failure  due  to  divergence  by  the 
nature  of  the  recorded  signal.  When  the  speed  was  increased  to  36  knots,  the  model  failed  in 
bending,  and  several  segments  broke  off,  as  shown  in  Figure  7.  The  strain  gage  output 
showed  no  oscillation,  establishing  that  the  failure  was  purely  static. 

It  was  expected  that  all  of  the  models  would  fail  similarly  at  about  36  knots,  unless  a 
flutter  condition  occurred  at  a  lower  «peed  as  for  Model  1.  Since  Model  3  was  quite  similar  to 
Model  2,  it  would  probably  not  flutter  below  36  knots  and  would  give  no  useful  information  by 
being  tested.  Model  4,  however,  had  a  substantially  lower  mass  ratio  than  Model  2  and  might 
possess  a  lower  flutter  speed.  Therefore,  Model  4  was  tested  up  to  30  knots,  with  the  result¬ 
ing  decay  factors  and  frequencies  shown  in  Figures  12  and  13.  As  in  the  case  of  Model  2,  the 
measured  decay  factors  were  very  high,  and  coul !  not  be  read  above  15  knots.  Model  excur¬ 
sions  due  to  flow  fluctuations  were  becoming  visually  noticeable  at  30  knots,  implying  close¬ 
ness  to  static  failure.  As  has  been  previously  described,  air  leakage  from  Model  4 
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necessitated  several  repairs  to  the  Silastic  seams.  The  value  of  further  testing  became 
more  doubtful  when  all  air  leakage  stopped  during  a  run  at  28  knots,  indicating  that  the 
model  had  filled  with  water.  A  60  Hz  signal  appeared  in  the  strain  gage  output  shortly  after¬ 
ward,  further  reducing  its  readability.  Since  the  water  had  changed  the  inertial  character¬ 
istics  of  the  model,  further  testing  was  abandoned  and  the  model  was  removed  from  the  tunnel 
for  retesting  of  its  structural  stiffnesses.  Subsequent  examination  of  the  model  showed  a 
number  of  openings  in  the  segments  along  the  spar  where  pieces  had  broken  off.  Apparently 
flooding  occurred  through  these  openings  in  the  segments  rather  than  through  the  Silastic 
seams. 


THEORETICAL  FLUTTER  ANALYSIS 

The  theoretical  analysis  was  intended  to  explore  the  validity  of  the  Yates  flutter 
theory  in  the  low  mass  range.  It  was  recognized  that  the  essential  feature  of  the  Yates 
theory  was  the  modification  that  was  made  to  two-dimensional  Theodorsen  loading  in  order 
to  represent  three-dimensional  flow  effects.  Accordingly,  the  three  calculations  described 
below  all  incorporate  spanwise  variations  in  loading. 

The  structural  representation  used  by  Ya.es  in  his  analysis  was  not  considered  to  be 
an  essential  feature  of  the  theory.  Therefore,  Yates’  original  six-mode  Rayleigh-type 
analysis  was  not  preserved  in  the  present  calculations.  Instead,  two  calculations  were  made 
at  NSRDC  using  a  lumped  parameter  structural  representation.  A  third  calculation  was  made 
at  SwRl  using  a  two-mode  Rayleigh-type  analysis. 

Both  calculations  made  at  NSRDC  were  corrected  for  structural  damping  as  follows. 

An  exponential  decay  factor  8S  based  on  the  measured  structural  damping  was  added  to  the 
calculated  decay  factor  8h.  The  factor  8s  was  given  by 

8s  =  C,  Ju2  +  Sh 

and  the  totai  decay  factor  by 

8  =  8h  +  8s 

No  structural  damping  was  included  in  the  calculation  made  at  SWRI. 

YATES  FLUTTER  THEORY 

The  first  theory  to  be  discussed  employed  the  original  Yates  loading  theory  and  will 
be  referred  to  as  the  Yates  flutter  theory  despite  a  different  structural  formulation  than  Yates 
used.  The  hydrodynamic  loading  of  Yates  consisted  of  stripwise  application  of  two- 
dimensional  Theodorsen  loading,  with  loading  due  to  circulatory  flow  modified  for  the  effects 


10 


of  finite  span.  The  modification  was  accomplished  by  substituting  three-dimensional  distri¬ 
butions  of  lift  slope  and  center  of  pressure  iocation  in  place  of  the  two-dimensional  values 
of  *2ff  and  quarter  chord,  respectively.  The  three-dimensional  spanwise  distributions  used 
were  the  same  values  used  bv  Vales  in  his  calculation  for  an  aspect  ratio  5  foil.  These 
values,  which  were  taken  from  Reference  11,  are  shown  in  Figure  14.  The  distributions 
were  calculated  with  lifting  surface  theory. 

The  modified  hydrodynamic  loading  was  applied  to  a  lumped  parameter  representation 
of  the  hydrofoil  structure.  The  SADSAM  HI  computer  program15,  !6,  17  was  used  to  generate 
structural  information  based  on  measured  model  parameters.  Each  model  was  represented  by 
structural  stiffness  and  inertial  parameters  at  10  equally  spaced  spanwise  stations  at  which 
hydrodynamic  loading  was  applied.  This  permitted  the  calculation  of  20  coupled  structural 
modes. 

Predicted  natural  frequencies  i..  vacuum  for  Models  1  through  4  and  the  parent  model 
are  given  in  Tables  2  through  6.  The  predicted  frequencies  range  from  5  percent  high  to  11 
percent  low,  averaging  5  percent  lower  than  experiment.  The  modes  are  labelled  in  Tables  2 
through  6  according  to  predominant  uncoupled  mode  components.  The  predicted  nodal  lines 
in  Figure  6  show  excellent  agreement  with  experiment  except  for  the  second  and  third  modes 
of  Model  4.  Th<s  discrepancy  may  have  been  caused  by  a  structural  irregularity  in  Model  4, 
as  is  suggested  by  the  nodal  line  for  the  second  mode.  It  was  concluded  that  the  structural 
representation  was  approximately  correct  for  Models  1,  2,  and  3,  but  not  for  Model  4.  Despite 
the  substantial  variation  in  model  densities,  it  is  interesting  to  note  that  the  mode  shapes  of 
all  the  models  are  nearly  the  same,  except  for  the  fifth  and  sixth  modes  of  Model  4  which  have 
interchanged  their  relative  positions. 

Using  the  zero-speed  or  noncirculatory  part  of  Yates’  hydrodynamic  loading,  the  nat¬ 
ural  frequencies  of  the  models  in  water  were  predicted  and  are  given  in  Tables  2  through  6. 

The  predicted  frequencies  range  from  1  percent  to  13  percent  below  experiment,  and  are  on 
the  average  9  percent  low.  This  result  implies  that  the  noncirculatory  hydrodynamic  loading 
is  not  well  represented  in  Yates’  treatment. 

Model  response  characteristics  in  the  form  of  exponential  decay  factors  were  calcu¬ 
lated  as  functions  of  speed  for  the  four  present  models  and  the  parent  model.  These  are  shown 
in  Figures  8  through  13  and  15  through  18  along  with  the  experimental  results.  Flutter  is 
predicted  to  occur  at  a  speed  at  which  the  exponential  decay  factor  becomes  zero.  The  pre¬ 
dicted  flutter  speeds  and  frequencies  are  indicated  in  each  figure. 

According  to  the  two  flutter  speeds  that  were  obtained  experimentally,  the  predictions 
of  the  Yates  theory  were  overconservative.  The  20.1-knot  prediction  for  Model  1  was  23  per¬ 
cent  below  the  observed  flutter  speed,  and  the  33.7-knot  prediction  for  the  parent  model  was 
30  percent  below  experiment.  Similarly,  the  26. 4-knot  flutter  prediction  for  Model  2  was  27 
percent  below  the  static  failure  speed  of  the  model.  The  magnitudes  of  the  predicted  decay 
factors  were  much  lower  than  the  observed  magnitudes  throughout  the  speed  range,  showing 
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that  damping  was  consistently  underestimated.  Frequency  predictions  were  consistently  be¬ 
low  experiment  at  low  speeds  and  above  experiment  at  high  speeds.  No  indication  was  evi¬ 
dent  in  the  calculations  of  the  decrease  in  flutter  speed  which  had  been  found  at  low  values 
of  mass  ratio  in  the  Yales  study. 

NSRDC  FLUTTER  THEORY 

Ti.e  second  theory  to  be  discussed  consisted  of  an  extension  of  the  Yates  theory  pre¬ 
viously  described.  The  extended  theory  modified  the  noncirculatory  loading  to  account  for 
three-dimensional  flow  effects.  This  extension  was  based  on  a  suggestion  made  by  Yates  in 
Reference  7. 

The  noncirculatory  loading  terms  were  multiplied  by  a  spanwise  function  which  was 
inferred  from  the  virtual  mass  expression  for  a  flat  plate.  This  modification  was  chosen  in 
view  of  the  presence  of  virtual  mass  and  moment  in  the  noncirculatory  loading  expression. 
The  function  consisted  of  the  fraction  of  two-dimensional  virtual  mass  of  the  foil  outboard 
of  the  spanwise  position  being  considered  and  was  given  by 

A* 

v/l  +  (A*)2 

In  this  expression  A*  is  the  local  aspect  ratio,  which  is  the  aspect  ratio  of  the  foil  outboard 
of  the  local  position,  when  a  reflecting  plane  exists  at  the  local  position.  The  distribution 
factor  p  is  shown  in  Figure  19  for  an  aspect  ratio  5  hydrofoil.  Total  virtual  mass  values 
from  this  function  agreed  well  with  measured  values  given  in  Reference  1. 

Predicted  natural  frequencies  in  still  water  for  the  five  models  are  given  in  Tables  2 
through  6.  The  predictions  range  from  10  percent  above  experiment  to  9  percent  below, 
averaging  1  percent  above  the  measured  values.  These  results  are  significantly  better  than 
the  predictions  of  the  unmodified  Yates  theory,  and  it  is  concluded  that  noncirculatory  load¬ 
ing  is  well  represented  by  the  spanwise  function  used. 

Model  flutter  characteristics  calculated  by  the  NSRDC  flutter  theory  are  shown  in  Fig¬ 
ures  8  through  13  and  15  through  18.  The  flutter  speed  predictions  of  the  NSRDC  flutter  theory  were 
conservative,  but  not  overcon servalive  as  were  those  of  the  Yates  flutter  theory.  The  pre¬ 
dicted  flutter  speed  of  23.1  knots  for  Model  1  was  7  percent  lower  than  experiment,  while  the 
Yates  theory  predicted  a  flutter  speed  which  was  23  percent  low.  Similarly,  the  41.8-knot 
prediction  for  the  parent  model  was  13  percent  below  experiment,  as  compared  to  a  Yates 
theory'  prediction  which  was  30  percent  low.  Flutter  speeds  for  the  remaining  models  were 
predicted  to  be  higher  according  to  the  NSRDC  theory  than  according  to  the  Yates  theory, 
although  the  prediction  for  Model  2  wrs  still  below  the  failure  speed  of  36  knots.  Frequency 
predictions  of  the  NSRDC  theory  were  fairly  accurate  over  the  speed  range  tested  and  showed 
slightly  better  agreement  with  experiment  than  the  Yates  theory  predictions. 
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In  addition  to  the  more  accurate  flutter  speed  predictions  of  the  NSRDC  theory,  the  de- 
cay  factors  were  more  accurately  predicted  over  the  entire  speed  range  in  comparison  with  the 
Yates  theory.  The  improvement  consisted  of  an  increase  in  the  predicted  decay  factors  to 
approximately  the  measured  values. 

The  predicted  decay  factors  and  flutter  speeds  of  both  the  NSRDC  theory  and  the  Yates 
theory  were  significantly  affected  by  inclusion  of  structural  damping.  This  raised  the  calcu¬ 
lated  flutter  speeds  from  10  to  35  percent  for  the  Yates  theory  and  from  6  to  14  percent  for 
the  NSRDC  theory.  It  is  therefore  concluded  that  structural  damping  was  not  negligible  for 
the  present  models  and  should  be  included  in  future  experimental  and  theoretical  studies. 

SWRI  FLUTTER  THEORY 

The  third  theory  to  be  discussed  is  another  modification  of  the  original  Yates  flutter 
theory,  and  is  described  as  Case  6c  in  Reference  11.  It  will  be  referred  to  as  the  SwRI 
flutter  theory.  In  this  calculation,  two-dimensional  Theordorsen  hydrodynamic  loading  was 
modified  to  incorporate  measured  spanwise  lift  slopes.  The  lift  slope  distribution  used  is 
shown  in  Figure  14a.  The  center  of  pressure  location  was  left  unchanged  from  its  quarter 
chord  position  corresponding  to  two-dimensional  flow.  The  hydrodynamic  loading  was  inte¬ 
grated  over  the  span  of  the  model  in  a  two-mode  Rayleigh-type  analysis.  No  structural  damp¬ 
ing  was  included. 

Flutter  characteristics  of  the  hydrofoil  models  calculated  by  the  SwRI  theory  are  shown 
in  Figures  8-13  and  15-18.  The  flutter  speed  predictions  agreed  very  closely  with  the  two 
experimental  flutter  speeds,  and  were  more  accurate  than  the  NSRDC  predictions.  Further¬ 
more,  the  prediction  of  48.0  knots  for  Model  2  was  above  the  speed  at  which  the  model  failed, 
so  that  the  SwRI  theory  may  not  be  as  conservative  at  low  mass  ratios  as  the  two  theories 
previously  discussed.  The  agreement  between  measured  and  calculated  decay  factors  was 
very  good.  Frequency  predictions  were  slightly  lower  than  the  NSRDC  theory  predictions. 

The  response  curves  for  the  first  mode  were  not  graphed,  but  were  found  to  be  stable. 

DISCUSSION 

The  flutter  speed  predictions  of  all  three  theories  have  been  plotted  in  Figure  20  as 
functions  of  mass  ratio.  All  theories  predicted  an  increasing  flutter  speed  as  mass  ratio 
decreased,  for  the  mass  ratios  of  the  models  tested.  The  flutter  curves  retained  their 
relative  positions  throughout  the  mass  ratio  range,  with  the  Yates  theory  predicting  the  low¬ 
est  flutter  speeds,  the  NSRDC  theory  predi  :ng  slightly  higher  flutter  speeds,  and  the  SwRI 
theory  predicting  the  highest  flutter  speeds. 

The  failure  of  the  three  flutter  theories  to  predict  a  flutter  speed  boundary  which  de¬ 
creased  as  mass  ratio  decreased  in  the  low  mass  ratio  range  was  unexpected  in  view  of  the 
Yates  prediction  in  his  original  analysis.7  The  decreasing  flutter  speed  found  by  Yates 


m  i  Hired  below  a  mass  rut io  of  approximately  0.3.  Differences  between  the  present  theories 
and  the  original  Vales  theory  may  have  accounted  for  the  discrepancy.  A  similar  downturn 
may  occur  for  the  present  theories  below  a  mass  ratio  of  0.202,  but  no  calculations  were  made 
to  determine  whether  tins  does  occur. 

Test  results  established  that  all  theories  were  conservative  at  mass  ratios  0.99  and 
0.9l!.‘t.  The  Yates  theory  and  the  NSRDC  theory  were  also  conservative  at  mass  ratio  0.42)5. 
The  other  predictions  could  not  be  evaluated  because  the  predicted  flutter  speeds  were  higher 
than  the  observed  bending  failure  speed. 

The  noncirculatory  loading  modification  used  in  the  NSRDC  flutter  theory  is  a  signifi¬ 
cant  improvement  over  the  unmodified  Yates  hydrodynamic  loading.  The  modification  pro- 
ductal  more  accurate  natural  frequency  predictions  for  the  hydrofoils  in  still  water  as  well  as 
improved  flutter  predictions.  The  effect  of  this  modification  on  flutter  predictions  at  very  low 
mass  ratio  was  not  determined. 

Although  the  SwRI  flutter  theory  gave  the  most  accuratc*  flutter  predictions,  it  has  less 
capability  than  the  NSRDC  flutter  theory  for  solving  flutter  problems  of  full-scale  hydrofoil 
systems.  The  deficiency  of  the  SwRI  theory  lies  in  its  use  of  only  two  structural  modes  of 
v  ibration.  Two  modes  cannot  adequately  represent  the  structural  properties  of  a  complex 
strut-foil  system.  The  theory  could  be  expanded  to  include  a  larger  number  of  modes,  as  was 
done  in  the  original  6-modo  analy  is  of  Yates.  However,  the  lumped  parameter  structural 
representation  used  in  the  NSRDC  theory  is  extremely  versatile  in  representing  complex 
structures,  including  pods  and  foils  attached  to  struts.  It  would  be  possible  to  combine  the 
lumped  parameter  representation  with  the  hydrodynamic  loading  used  in  the  SwRI  theory,  but 
this  should  await  further  experimental  confirmation  of  the  SwRI  theory. 

THEORETICAL  DIVERGENCE  ANALYSIS 

The  structural  failure  of  Model  2  at  36  knots  indicated  a  need  for  further  analysis  of 
the  divergence  characteristics  of  the  flutter  models.  Model  design  had  been  based  on  two- 
dimensional  load  theory,  which  had  predicted  .  center  of  pressure  at  the  quarter  chord  of  the 
model  and  consequently,  an  infinite  divergence  speed.  Although  a  divergence  failure  had  not 
been  predicted,  structural  failure  was  predicted  to  occur  as  a  result  of  bending  stress  as  a 
function  of  foil  angle  and  speed.  For  example,  a  uniform  foil  angle  of  3.6  deg  or  a  tip  twist 
angle  of  6.8  deg  would  have  caused  the  spar  to  yield  at  36  knots.9  Since  model  trim  and  tip 
excitations  had  been  held  to  much  smaller  values  than  these,  the  observed  failure  cannot  be 
explained  by  this  prediction.  It  is  concluded  that  the  two-dimensional  load  theory  is  inade¬ 
quate  to  predict  the  loading  on  the  foil. 

A  more  accurate  prediction  of  failure  characteristics  could  use  three-dimensional  lift- 
ing  surface  theory,  which  predicts  a  center  of  pressure  location  forward  of  the  quarter  chord 
on  finite  aspect  ratio  foils.  A  forward  location  of  the  center  of  pressure  would  increase  the 
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probability  of  either  divergence  or  bending  failure  by  causing  torsional  deflections  to  increase, 
thereby  increasing  the  stress  on  the  model. 

Therefore,  the  Yates  hydrodynamic  loading  formulation,  previously  used  in  the  Yates 
flutter  theory  and  based  on  a  three-dimensional  lifting-surface  calculation  of  lift  slope  and 
center  of  pressure,  was  used  to  calculate  model  failure  characteristics.  The  steady  loading 
terms  were  combined  with  the  lumped  parameter  stiffness  matrix  used  in  the  flutter  calcula¬ 
tion.  In  matrix  form,  resulting  equations  for  the  10  span.vise  stations  used  are 

[Kh]  =  IP0  {«,-)! 

i  =  1,  2,  ...  10  (1) 

[  Ke]  =  IV«,)i 

where  I/f^]  =  bending  stiffness  matrix, 

[ Kq]  =  torsional  stiffness  matrix, 

P0(a,-)  =  steady  lift  on  foil  due  to  angle  of  attack, 

M0(a,)  =  steady  moment  on  foil  due  to  angle  of  attack, 

A,.  =  local  translational  displacement  of  foil  elastic  axis  from  unstressed  position, 

=  local  torsional  displacement  of  foil  elastic  axis  from  unstressed  position,  and 
a£  =  local  angle  of  attack  between  foil  and  free  stream. 

The  coordinates  A  and  0  correspond  to  elastic  displacements  of  the  foil  spar.  They  a<-e  re¬ 
ferred  to  the  unstressed  position  of  the  foil,  not  to  the  flow  direction,  and  are  zero  when,  for 
example,  the  foil  is  motionless  in  ctill  water.  On  the  other  hand,  the  angles  <*,-  are  defined 
as  the  angles  that  local  chords  make  with  the  free  stream.  Thus  the  a,-  include  both  struc¬ 
tural  twist  angles  and  misalignment  angles  between  the  foi!  and  the  flow  direction.  The  com¬ 
posite  nature  of  the  a,-  can  be  expressed  by  the  relation 

ai  =  0;  +  Y- 

where  y,-  =  local  angle  of  unstressed  foil  relative  to  the  free  stream.  Relationships  among  the 
three  angular  coordinates  are  shown  in  Figure  21. 

Since  the  lift  and  moment  expressions  are  proportional  to  angle  of  attack  a £,  Equation 
(1)  may  be  written 


[Kh]  IA,.l  =  \P0(6.)\  +  IPq  (y,)l 

I**]  10,-1  =  \M0(6.)\  +  \M  0  (>',-)! 

The  terms  in  the  second  of  the  above  equations  can  be  combined  to  give 

[Ke-M]  |0,-i  =  Uf0(y,.)!  (2) 
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where*  l*/  i  is  a  diagonal  matrix  defined  by 

u7 1  i#  i  =  i .</„(<?,.) i 

A  non-lmial  solution  to  the  homogeneous  form  of  Equation  (2)  exists  when  the  determinant  of 
the  coefficient  matrix  vanishes.  The  determinant  will  vanish  for  some  value  or  values  of  free 
stream  velocity  V /t.  U is  the  divergence  speed  because  the  vanishing  of  the  determinant  of 
[A'^  -  1/1  implies  that  all  0.  are  infinite  for  a  non-zero  set  of  values  I M ( y,) I.  For  speeds  be¬ 
low  l  1 Equation  (1)  can  be  solved  for  specified  initial  flow  misalignment  angles  v.  to  give 
the  bending  and  torsional  displacements  of  the  foil.  Stress  values  can  then  be  determined 
from  the  displacements. 

The  divergence  speed  of  Model  2  was  calculated  using  the  stiffness  matrix  [A'^l  and 
the  steady  part  of  the  hydrodynamic  loading  which  had  been  previously  used  in  the  Yates  and 
NSRDC  flutter  calculations.  A  value  of  44.0  knots  was  obtained  for  V q.  Since  this  speed 
was  substantially  higher  than  the  observed  failure  speed,  further  calculations  were  made  to 
determine  whether  the  model  might  have  failed  in  bending  prior  to  reaching  divergence. 

The  response  of  a  typical  hydrofoil  to  hydrodynamic  loading  near  the  divergence  speed 
has  been  discussed  in  Reference  18.  The  response  of  the  foil  varies  from  zero  at  zero  speed 
to  infinity  at  the  divergence  speed.  Foil  torsional  deflection  may  be  described  in  terms  of  a 
magnification  effect  in  which  ratio  of  the  model  tip  twist  angle  to  the  misalignment  angle 
varies  from  zero  to  infinity.  When  the  ratio  is  unity,  a  given  misalignment  angle  produces  an 
equal  tip  twist.  This  magnification  ratio  was  calculated  using  the  lumped  parameter  repre¬ 
sentation  for  Model  2  and  is  shown  in  Figure  22.  The  misalignment  angle  was  assumed  to  be 
constant  along  the  span  of  the  model  for  purposes  of  calculation.  Unit  magnification  of  misa¬ 
lignment.  ang!  ?  occurs  at  26.5  knots,  or  60  percent  of  the  divergence  speed.  At  36  knots  it  is 
predicted  that  the  structural  twist  of  the  model  tip  will  be  3.5  times  the  misalignment  angle. 

Correspondingly  high  bending  displacements  are  predicted  to  accompany  the  torsional 
displacements.  In  order  to  determine  the  approximate  misalignment  angle  that  would  cause 
failure  in  bending,  a  spar  yield  boundary  was  calculated  using  the  yield  stress  of  the  spar, 
350,000  lb/ sq  in.  This  boundary  is  plotted  in  Figure  23.  A  misalignment  angle  of  0.92  deg 
was  predicted  to  cause  bending  failure  at  36  knots. 

Motion  pictures  of  Model  2  showed  that  ambient  flow  fluctuations  caused  tip  deflections 
of  as  much  as  1  deg  in  torsion  and  0.4  in.  in  bending,  about  the  mean  position,  at  35  knots. 
The  flow  fluctuation  angles  were  not  measured  directly. 

On  the  basis  of  calculated  model  response  and  observed  model  deflections  it  is  con¬ 
cluded  that  Model  2  failed  structurally  in  bending  at  36  knots.  The  failure  was  caused  by  too 
close  an  approach  to  the  divergence  speed  of  the  model.  The  failure  was  initiated  by  small 
flow  fluctuations  in  the  tunnel  jet  which  were  magnified  by  the  model  divergence  character¬ 
istics.  Each  of  the  other  models  would  have  failed  similarly  at  nearly  the  same  speed  be¬ 
cause  of  their  similar  elastic  properties. 
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Tilt'  divergence  theory  used  agrees  reasonably  well  with  observed  model  behavior. 

Future  model  design  should  be  based  on  a  similar  divergence  calculation.  In  order  to  test  a 
model  successfully  at  a  given  speed  the  divergence  speed  must  be  substantially  higher  than 
the  test  speed.  An  appropriate  test  speed  limit  might  be  60  percent  of  the  calculated  diver¬ 
gence  speed,  at  which  speed  the  response  of  Model  2  equalled  the  flow  fluctuation  angle.  At 
this  speed  model  response  to  either  flow  fluctuations  or  deliberate  excitation  would  be  limited 
to  the  initial  disturbance  amplitude  and  would  cause  little  interference  with  damping  measure¬ 
ments.  Additional  calculations  would  be  required  to  ensure  that  the  model  would  not  fail  in 
bending  as  a  result  of  the  model  excitation. 

The  divergence  speed  of  models  of  the  present  type  could  be  increased  by  placing  the 
elastic  axis  farther  forward,  in  order  to  have  the  center  of  pressure  of  the  loading  aft  of  the 
elastic  axis.  Figure  24  illustrates  the  dependence  of  theoretical  divergence  speed  on  elastic 
axis  position  for  Model  2.  If  the  elastic  axis  were  located  at  a  =  -0.58,  corresponding  to  an 
elastic  axis  at  21  percent  of  the  chord  aft  of  the  leading  edge,  the  predicted  divergence  speed 
would  be  83.3  knots  and  a  test  speed  of  50  knots  would  equal  60  percent  of  the  divergence 
speed. 

CONCLUSIONS 

Four  hydrofoil  models  were  constructed  which  successfully  scaled  a  previously-tested 
hydrofoil  to  smaller  size  and  to  lower  mass  ratios.  Flutter  was  obtained  for  the  model  of 
mass  ratio  0.963.  The  models  of  mass  ratio  0.455,  0.395,  and  0.202  were  subject  to  static 
failure  in  bending  before  flutter  was  encountered. 

The  results  for  Model  1  demonstrate  that  it  is  possible  to  approach  a  hydrofoil  flutter 
condition  very  closely  without  destroying  the  model.  It  is  concluded  that  the  flutter  inception 
speed  of  a  hydrofoil  can  be  determined  in  some  cases  without  the  occurrence  of  flutter  if  damp¬ 
ing  measurements  are  made  while  speed  is  increased,  and  if  speed  increments  are  kept  small. 

The  flutter  characteristics  of  the  models  were  calculated  using  three  versions  of  the 
Yates  hydrodynamic  loading  theory.  Both  the  Yates  flutter  theory  and  the  NSRDC  flutter 
theory  showed  a  tendency  to  be  overconservative.  Inclusion  cf  structural  damping  in  these 
calculations  had  a  significant  effect  on  calculated  flutter  speeds.  The  SwRI  flutter  theory 
agreed  well  with  experiment  at  mass  ratios  of  0.933  and  0.99,  but  could  not  be  confirmed  at 
lower  values  of  mass  ratio. 

A  spanwise  modification  to  the  noncirculatory  loading,  contained  in  the  NSRDC  theory, 
significantly  improved  the  predictions  of  natural  frequencies  in  water  and  of  model  damping 
characteristics  and  flutter  speeds.  Additional  calculations  are  required  to  determine  the 
nature  of  the  low  mass  ratio  flutter  boundary  with  and  without  this  modification. 

The  static  bending  failure  observed  on  the  mass  ratio  0.455  model  was  caused  by  too 
close  an  approach  to  the  divergence  speed  of  the  model.  The  failure  was  initiated  by  small 
flow  fluctuations  in  the  tunnel  jet  which  were  magnified  by  the  model’s  divergence  character¬ 
istics.  Each  of  the  other  models  would  have  failed  simiiarly  at  nearly  the  same  speed  because 
of  th'»  similar  elastic  properties  of  the  models.  The  observed  failure  was  consistent  with  the 
calculated  divergence  speed  of  44  knots. 
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RECOMMENDATIONS 


1.  The  present  tests  were  unsuccessful  in  defining  flutter  boundaries  at  low  mass  ratio 
because  of  structural  failure  of  the  models.  Future  small-scale  flutter-inodel  design  should 
employ  three-dimensional  hvdrodynamic  load  theory  in  calculating  the  spanwise  center  of 
pressure  location  for  divergence  speed  predictions.  Test  speeds  should  be  limited  to  60  per¬ 
cent  of  the  divergence  speed  to  eliminate  undue  model  deflections  and  to  preclude  st-uctural 
failure. 

‘2.  Full-scale  hydrofoil  craft  utilize  low  mass  ratio  struts  with  relatively  massive  pods 
and  foils  attached  to  the  ends  of  the  struts.  It  has  been  shown  that  flutter  speeds  of  struts 
can  be  radically  lowered  by  such  tip  elements.1,  19  The  effects  of  foil  configuration,  pod 
inertia  and  shape,  and  submergence  depth  should  be  further  investigated  both  experimentally 
and  theoretically.  Addition  of  tip  masses  would  permit  flutter  testing  of  models  similar  to 
those  described  in  this  report. 

.‘5.  Future  flutter  testing  should  utilize  damping  measurement  to  detect  the  approach  of 
flutter.  Methods  of  damping  measurement  should  be  developed  for  eventual  use  on  full-scale 
hydrofoil  craft. 
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Figure  3  —  Flutter  Model  and  Model  Support  Enclosure  Mounted  in  36-lnch  Variable  Pressure  Water  Tunnel 
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Figure  4  -  Section  through  Turntable 
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Figure  6  -  Thooietical  and  Experimental  Nodal  Lines  for  the  First 
Six  Natural  Frequencies  of  Flutter  Models  in  Air 
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MODE  1 
EXP:  15.7  Hz 
TH:  15.3  Hi 


TH:  31  83  Hi  (™1  TH:  55.68  Hi  (-) 


TH:  88  89  H  t  ( — I 


TH:  131.6  Ha  { — ) 


Figure  6a  —  Model  1 


EXP:  46.4  Hi  (•) 
TH:  43.45  Hi  (—) 


EXP:  79.8Ki(#) 
TH:  77.03  Hi  (-) 


EXP:  128.5  Hi 
TH:  122.1  Hi  (-) 


EXP:  181.3  Hi 
TH:  180.0  Hi  (— ) 


MODE  6 
EXP:  157.1  Hi 
TH:  140.2  Hi  (—) 


MODE  6 
EXP:  191.7  Hi 
TH:  194.2  Hx  ( — ) 


Figure  6b  —  Model  2 


E»P:  17.2  Hi  EXP:  49.2  Hi  (#)  EXP:  84.8  H.  EXP:  '42.0  H.  EXP:  200.3  H:  EXP:  214.5  Hi 

TH:  14.45  Hi  TH:  44.69  Hi  <—)  TH:82  45Hi(— )  TH:  130.4  Hi  (—)  TH:  192  2  H- { — )  TH:  200.2  Hi  (—) 

Figure  6c  —  Model  J 


MODE  1  MODE  2  MODE  3  MODE  4  MODE  5  MODE  6 

EXP:  22.5  Hi  EXP:  68.9  Hi  (#)  EXP:  105.9  Hi  (#)  EXP:  173  9  Hi  EXP: -  EXP: - 

TH:  23.04  Hi  TH:  61  02  Hi  (— )  TH:  108.4  Hi  (—)  TH:  179  4  Hi  ( — )  TH:  261  6  Hi  (— )  TH:  268  3  Hi  (-) 


Figure  6d  —  Model  4 


Figure  7  -  Model  2  after  Static  Failure  in  Bending 
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Figure  8  —  Exponential  Decay  Factor  S  as  a  Figure  9  —  Frequency  as  a  Function  of 

Function  of  Speed  for  Model  1  Speed  for  Model  1 
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Figure  14  —  Spanwise  Load  Distributions  for 
Aspect  Ratio  5  Rectangular  Foil 


FREQUENCY  IN  HERTZ 
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Figure  16  -  Frequency  as  a  Function  of  Speed  for  Model  3 


gure  17  —  Exponential  Decay  Factor  5  as  a  Function  Figure  18  —  Frequency  as  a  Function 

of  Speed  for  the  SwRI  30-Inch  Flutter  Model  of  Speed  for  the  SwRI  30-Inch 

Flutter  Model 


NONCIRQJLATORY 


Figure  19  -  Noncirculatory  Modification  Factor  p  as  a  Function  of 
Spanwise  Position  for  an  Aspect  Ratio  5  Hydrofoil 


Figure  20  —  Theoretical  and  Experimental  Flutter  Speeds  for 
Models  1-4  as  a  Function  of  Mass  Ratio 


Figure  21  -  Foil  Deflectior  Angles 
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igure  23  —  Misalignment  Angle  y  to  Cause  Spar  Failure 
in  Bending  for  Model  2  as  a  Function  of  Speed 
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0  4  0  2  0  -0  2  -0  4  -0  6 

ELASTIC  AXIS  POSITION. 

Figure  24  -  Theoretical  Divergence  Speed  as  a  Function  of 
Elastic  Axis  Position  for  Model  2 
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TABLE  l 


Structural  Parameter  Values  for  the  Flutter  Models 


SWRI  30-Inch 

Parent  Made! 

Model 

1 

gg 

Model 

4 

0.95 

0.963 

0.455 

0.395 

0.202 

0.524 

0.524 

0.523 

0.528 

0.523 

0.512 

0.508 

0.503 

0.511 

0.506 

El  (lb-in.2) 

3.40  x  106 

56,900 

55,300 

55,800 

56,700 

GJ  (lb-in.2) 

0.973  x  106 

19,200 

16,420 

16,440 

15,630 

/.  (in.) 

30 

15 

15 

15 

15 

TABLE  2 


Theoretical  and  Experimental  Natural  Frequencies  for  SwRI  30-Inch  Parent  Mode! 

H  =  0.90 


Made  Number 
and  Type 

B  =  Bending 

T  -  Torsion 

Theoretical 
(in  Vacuum) 

Experimental 

(in  Air) 

Theoretical 
(in  Water) 

Experimental 
(in  Water) 

Yates 

Hydrodynamics 

NSRDC 

Hydrodynamics 

1  (IB,  IT) 

10.07 

9.63 

7.17 

8.03 

7.3 

2  (IT,  IB) 

28.28 

31.9 

22.15 

23.57 

23.8 

3  (2T,  IT,  2B,1B) 

50.17 

- 

37.51 

39.96 

- 

4  (3T,  2T,  2B) 

79.35 

- 

59.05 

62.29 

- 

5  (4T,  3T,  2T,  2B) 

116.7 

- 

87.90 

92.18 

- 

6  (5T,  2T,  2B) 

126.4 

- 

97.75 

103.1 

- 

j  NOTE  All  frequencies  are  given  in  Hertz. 
I - 
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TABLE  3 


Theoretical  and  Experimental  Natural  Frequencies  for  Mode?  1 

It  =  0.963 


: 


S 


Mode  Number 
ond  Type 

B  Bending 

7  Torsion 

Theoretical 
(in  Vocuum) 

Experimental 
(in  Air) 

Theoretical 
(in  Water) 

Experimental 
(in  Woter 

Electro- 

mognetic 

Shoker 

.  Tip 
Excitation 

Yates 

Hydro- 

d/nomics 

NSRDC 

Hydro¬ 

dynamics 

Tip 

Excitation 

1  (IB,  IT) 

10.78 

11.1 

11.0 

7.59 

8.19 

8.4 

2  (IT,  IB) 

31.83 

35.1 

33 

24.76 

24.12 

27 

3  (2T,  IT,  23,  iB) 

55.68 

58-61 

- 

•  51 

41.12 

- 

4  (3T,  2T,2B) 

88.89 

99.0 

- 

65.58 

68.38 

- 

(4T,  3T,  2T,  2 B)  Vacuum 

5  (2T,  2B)  Woter 

131.6 

144.8 

- 

98.54 

98.50 

- 

(5T,  2T,  2B)  Vocuum 

6  (4T,4B,  2B)  Woter 

140.2 

157.1 

- 

107.6 

102.9 

--  - 

NOTE:  All  frequencies  are  g.ven  in  Hertz. 


TABLE  4 


Theoretical  and  Experimental  Natural  Frequencies  for  Model  2 

It  =  0.455 


Mode  Number 
ond  Type 

B  =  Bending 

T  =  Torsion 

Theoretical 
(in  Vacuum) 

Experimentol 
(in  Air) 

Theoretical 
(in  Water) 

Experimental 
(in  Water) 

Tip 

Excitation 

Yates 

Hydro- 

dynomics 

NSRDC 

Hydro¬ 

dynamics 

Tip 

Excitation 

1  (IB,  IT) 

15.30 

15.7 

15.6 

8.63 

10.13 

9.9 

2  (IT,  IB) 

43.65 

46.4 

42 

28.17 

31.10 

30 

3  (2T,  IT,  2B,  IB) 

77.03 

79.8 

- 

46.59 

51.68 

- 

4(3T,  2T,  2B) 

122.1 

128.5 

- 

73.20 

79.38 

- 

5  (4T,  3T,  2T,  2B) 

1S0.0 

181.3 

- 

109.7 

117.9 

- 

6(5T,2T,2B) 

194.2 

191.7 

.  " _ 

123.5 

134.5 

- 

NOTE:  All  frequencies  are  given  in  Hertz, 

Table*  r. 


I'hiMiri'i u-al  ami  Experimental  Natural  Frequencies  for  Model  •! 

p  -  D..J95 


Mod#  Number 
and  T  ype 

/•  Bending 
/  fcrsion 

T  h^oreti  cal 

m  Vacuum1 

Experimental 
'•  in  Aifl 

Theoretical 
(in  Water) 

Experimental 
(in  Water) 

Electro¬ 

magnetic 

Shaker 

-  T*P 

Excitotion 

Y  ates 
Hydro¬ 
dynamics 

1 

.  T'P 
Exctation 

i  IB,  IT 

16.45 

17.2 

- 

8.84 

10.55 

- 

2  IT  IB 

46.69 

49.2 

- 

28.95 

32.37 

- 

3  2T  IT,  2B,  IB 

82.45 

86.8 

- 

47.81 

53.32 

- 

4  3T.2T.2B' 

13C.4 

142.0 

_ 

74.97 

81.90 

- 

5  4T.3T.  2T,  2B 

192.2 

200.3 

- 

112.3 

121.3 

- 

6  5T,2T,  2B1 

208.2 

214.5 

_ 

- 

127  0 

140.5 

- 

NOTH  All  I'requ 

L 

eru  t **s  iirt*  given  in  Hert7. 

TABLE  6 

Theoretical  and  Experimental  Natural  Frequencies  for  Model  4 

H  =  0.202 


Mode  Number 
and  Type 

H  Bending 
/  Tarsian 

Theoretical 
(in  Vacuum) 

~t 

Experimental 
(in  Air) 

Theoretical 
(in  Water) 

Experimental 
(in  Water) 

Electro¬ 

magnetic 

Shoker 

Tip 

Excitation 

Yates 

Hydro¬ 

dynamics 

NSRDC 

Hydro¬ 

dynamics 

Tip 

Excitation 

1  IB,  IT 

23.04 

22.5 

- 

9.57 

11.77 

- 

2  IT,  18! 

61.02 

68.9 

71 

31.54 

35.90 

35 

!  3  2T,  IT  2B,  IB) 

108.4 

105.9 

- 

51.89 

59.67 

- 

j  4f3T.2T.2B 

i 

179.4 

173.9 

- 

80.62 

89.69 

- 

1  ■  4 T ,  2 B 1  Vacuum 

5  4T,3T,  2T,  2B)  Water 

261.6 

- 

- 

120.6 

132.6 

- 

3T,  2T,  IT,  2Bi  Vacuum 
!  6  (5T,2T,28)  Water 

268.3 

- 

- 

139.0 

155.9 

- 

NOTK  Al!  frequencies  «»rr  given  in  Her 

7. 

:ih 


TABLE  7 
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